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ABSTRACT 

Wc show model light curves of superluminous supernova 2006gy on the assumption 
that the supernova is powered by the collision of supernova ejecta and its dense cir- 
cumstellar medium. The initial conditions are constructed based on the shock break- 
out condition, assuming that the circumstellar medium is dense enough to cause the 
shock breakout within it. We perform a set of numerical light curve calculations by 
using a one-dimensional multigroup radiation hydrodynamics code STELLA. We suc- 
ceeded in reproducing the overall features of the early light curve of SN 2006gy with 
the circumstellar medium whose mass is about 15 Mq (the average mass-loss rate 
~ 0.1 Mq yr~^). Thus, the progenitor of SN 2006gy is likely a very massive star. The 
density profile of the circumstellar medium is not well constrained by the light curve 
modeling only, but our modeling disfavors the circumstellar medium formed by steady 
mass loss. The ejecta mass is estimated to be comparable to or less than 15 Mq and 
the explosion energy is expected to be more than 4 x 10^^ erg. No ^^Ni is required 
to explain the early light curve. We find that the multidimensional effect, e.g., the 
Rayleigh- Taylor instability, which is expected to take place in the cool dense shell 
between the supernova ejecta and the dense circumstellar medium, is important in 
understanding supernovae powered by the shock interaction. We also show the evolu- 
tion of the optical and near-infrared model light curves of high-redshift superluminous 
supernovae. They can be potentially used to identify SN 2006gy-like superlmninous 
supernovae in the future optical and near-infrared transient surveys. 

Key words: supernovae: individual (SN 2006gy) — circumstellar matter — stars: 
mass-loss — early Universe 



1 INTRODUCTION 



Recent untargctcd transient surveys such as Texas Super- 
nova Sear ch (Quimby 2006a), Palomar Transient Factory 
l|Law et al.i2009l:lRau et al.| 2009:'). Catarina Real-time Tran- 
sient Survey (|Drake et al. " ]|200!D . Panoramic Survey Tele- 



scope & Rapid Response System (|Hodapp et al.ll2004 ). dis- 
covered new kinds of supernovae (SNe). Among the most 
specutaclar are superluminous SNe (SLSNe) which be- 
come typica l ly br ig hter than ~ — 21 m a g in optical (e.g 



Quimby et al.l 



— 21 m ag 
] l2007l, B 



20 111: iRest et al 



Smi th et all |2007| : 
2011: 'Miller et al.' '2009; 'Gezari et al."2009': 'Gal- Ya m et al 
12009. : .Drake et al.. .20^ . ,2011, : ,Chatzopoulos et al.il2oTlI : 
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Pastorello et al. 201C 


: Barbarv et alj 20091: Chomiuk et alj 


201ll:lLeloudas et all 


2OI2I: iBereer et al.ll2012l). 



SN 2006gy is the first example of SLSNe with detailed 
photomentric and spectroscopic observations. The luminos- 
ity of SN 2006gy stayed brighter than —21 mag in the R 
band for more than 50 days (see Section [2] for a brief sum- 
mary of observations). It was the most luminous SN ever 
reported at that time. Many theoretical models have been 
proposed to explain the huge luminosity of SN 2006gy. SN 
2006gy was, at first, hnked to Type la SNe exploding in- 
side H-rich circums t ellar medium (CS M ) (e.g., SN 2002ic; 



iHamuv et aO l2003l : IPeng et al.l l2004l l jofek et al.1 l200'if ). 

However, the total observed radiation energy exceeded the 
available energy from a Type la explosion and thus the no- 
tion was adandoned. 

Currently, SN 2006gy is thought to be caused by the 
death of a massive star. Several models have been proposed 
to explain the observed features as follows. All of them posit 
a massive progenitor star. 



rather than an SN explosion in this case. This kind of lumi- 
nous transients are called pulsational pair-instability SNe. 

(iv) Spin-down of newly born magnetars 

If a magnetar born at the time of the core collapse of 
a massive star has suitable magnetic field and spin, it can 
release its rotational energy efficiently just after the SN ex- 
plosion. If the radiation energy released at the magnetar 
can be converted to thermal energy, SN ejecta surround- 
ing the ma gnetar is heated up and the SN can become 
very brigh t (Kaseii fc Bildsten|[2oTol : IWooslevllioiol . see also 
iMaeda et al. 2007). 

(v) Transition of a neutron star to a quark star 
Neutron in a neutron star formed at the core collapse 

can be further discomposed to quarks and they can form 
a quark star. When a neutron transforms to quarks, la- 
tent heat due to the phase transition can heat up the SN 
ejecta surrounding the quark star and make an SLSN (e.g. , 
iLeahv fc Ou"vedll2008l:IOu ved et al.l l20ll : iKostka et al1l2012l . 
see also iBenvenuto fc Lu goncs 199^). 



(i) Large production of ^''Ni 

The energy released by the radioactive decay of ^®Ni is 
a common luminosity source of SNe. However, more than 
10 Mq of ^®Ni would be req uired to account for t he peak lu- 
minosity of SN 2006gy fe.g.. lNomoto et al.llioOTi '). The large 
amount of ^^Ni is not easil y produced even by ver y high 
energetic core-collapse SNe ( Umeda fc Nomotoll2008l'). Pair- 



instability SNe (PISNe, e.g.. iBarkat. Rakavv. fc Sack|[l967l : 



[Rakavy fc Shaviv 1967) can produce the large ^^Ni but the 
rising time of SN 2006gy is foun d to be much shorter than 
that expected from a PISN (e.g., Kasen. Wooslev. fc HegeJ 
I2OIII : iNomoto et al] l2007l : iMoriva et al.ll2010l V The weak- 
ness of Fe lines in spectra at the late phase is a lso against 
the ^®Ni heating scenario |Kawabata et al. 1 |2009l ). 

(ii) Interaction between SN ejecta and dense CSM 
The interaction between SN ejecta and dense CSM cre- 
ated by the progenitor before the explosion can make the 
luminosity of an SN very large. Strong shocks convert the 
kinetic energy of SN ejecta to ther mal energy which i s even - 
tually released as radiation energy. ISmith fc McCravl (|2007l ) 
suggest that the declining phase of the light curve (LC) of 
SN 2006gy can be e xplained by emission f rom a shocked 
~ 10 Mq CSM shell. lAgnoletto et al.l l|2009l) combined the 
^^Ni and shock interaction to explain the early LC of SN 
2006gy. They argue that the rising part of LC is powered 
by ®^Ni decay and the shock interaction comes into play 
from around the LC peak. Then, the required ^®Ni mass is 
reduced to a few Mq. 

(iii) Pulsational pair-inst a bility 

IWooslev. Blinnikov. fc HegeJ |2003) relate SN 2006gy to 
the pulsational pair-instability. A very massive star with 
mass ~ 80 — 130 Mq at the zer o-age main sequence un- 
dergoes strong pulsa tions (see also lUmeda fc Nomotdliool : 
lOhkubo et al.ll2009l ). which induce the intermittent exten- 
sive mass loss. If the mass ejected by such an eruption is 
caught up by the mass released by the next eruption, the two 
massive shells collide and the star can eventually be as lumi- 
nous as SN 2006gy. In this case, the luminosity source is the 
kinetic energy of the secondly ejected materials. Hence, the 
radiation mechanism is essentially the same as that of (ii) the 
interaction between SN ejecta and dense CSM. The ejecta 
from inside is released by the pulsational pair-instability, 



Among the possible heating sources, the CSM interac- 
tion (ii and iii) is plausible because the spectra of SN 2006gy 
show narrow P-Cygni profiles that indicate the existence of 
a CSM outflow with ~ 100 km s^'^. Although the observed 
low X-ray flux of SN 2006gy (see Section [2]) might appear to 
contradict the interaction scenario, the CSM could be op- 
tically thick enough to absorb X-rays from the shock wave , 
or the shock wa ve itself could be cool (e.g.. Blinnikovl l2008l : 
IChevalier fclrw in 2012; Svirski, Nakar, fc Sari 2012). 

Although the scenario (ii), i.e., the shock inter- 
action between SN ejecta and dense CSM, is sug- 
gested to be able to explain the early LC of SN 
200 6gy, detailed LC models have not been explore d 



yet. Chevalier fc IrwinI (I2OIII ): iMoriva fc Tominagal (|2012h : 
I Svirski. Nakar. fc Saril (|2012l ) investigate the possibility of 
explaining SN 2006gy by the interaction of SN ejecta and 
dense CSM alone. They assume that the CSM is opti- 
cally thick enough to cause the shock breakou t within it 



IWeaveilll976l: iNakar fc Sarill2010l . see also lOfek etal] 
iBalberg fc LoebT 201 ll about the shock break out in 
dense CSM). IChatzopoulos. Wheeler, fc Vinkd l|2012l ') show 
the bolometric LC from the interaction scenario based on 
their semi-analytic LC model (see Section [6.4p . The previ- 
ous works regarding the collision of SN ejecta and dense 
CSM so far analyse the LC of SN 2006gy based on analytic 
approaches with several simplifications. The detailed numer- 
ical modeling of the LC of SN 2006gy has not been done yet 
an d it is clearly required (see also the very recent attempt 
bv lGinzburg fc Balberej|2012l ). 

In this paper, we study model LCs for SN 2006gy by 
using numerical radiation hydrodynamics code STELLA. Our 
purpose of the modeling is to show numerically that the 
simple analytic shock breakout model actually works well to 
explain overall features of the early LC of SLSN 2006gy. We 
focus on SN 2006gy because it is the best observed SLSNe 
currently reported. We do not aim at the perfect fitting 
of the LC because it does not necessarily lead us to the 
correct physical parameters of SN 2006gy. This is because 
STELLA makes several simplifications to numerically treat 
the radiation hydrodynamics. We rather concentrate on the 
overall features of the LC of SN 2006gy and show that the 
dense CSM configurations predicted by the shock breakout 
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model are actually working well to reconstruct the LC of SN 
2006gy. 

The rest of the paper is organized as follows. Section 
[2] is a short summary of the observations of SN 2006gy. In 
Section [3l STELLA, which is used for our LC modeling, is 
summarized briefly. How multidimensional effect is included 
in one-dimensional code STELLA is also overviewed in Sec- 
tion [S] Our models are introduced in Section [4] and the LC 
calculations based on the models are shown in Section [S] 
Discussion is given in Section [B] We close the paper with 
conclusions in Section [T] We use the standard cosmology 
with Ho = 70 km s"^ Mpc~\ Q.m = 0.3, and f^A = 0.7 
when it is required. 



2 BRIEF SUMMARY OF OBSERVATIONS 

We briefly summarize the observational properties of SN 
2006gy. SN 2006gy was discovered by Texas Supernova 
Search on September 18, 2006 (U T) near the nu cleus of 
an early- type galaxy NGC 1260 (|Quimbvl [2"006bl '). There 
are sev eral suggested values for the extinction by the host 
galax y (|Ofek et al.ll2007l : ISmith et al.ll2007l : lAgnoletto et all 
20091). We adop t E{B - y)hoat = 0.40 mag following 
Agnoletto et al.l (|2009l ) wit h the Milky Way extinction 
E{B — V)mw = 0.16 mag (|SchIegel. Finkbeiner. fc Davij 
Il998l 'l. Thus, the total extinctio n \s E(B-V) = 0.56 mag or 



Ar = 1.3 mag with Rv = 3.1 (ICardelU. Clayton. &: Mathij 
Il98g| ). The distance modulus fi o f the host galaxy is also 
taken from lAgnoletto eFal] (|2009l ') (fj. = 34.53 mag). 

Follow up spectral observations classified SN 2006gy as 
Type Tin because o f the narrow H emission lines p resented 
in the spectra fsee ISchlegellll990l : iFil ippcnko T997I . for the 
details of Type Iln). The luminosity of SN 2006gy kept rising 
until October 25, 2006 (UT) and the peak R band luminos- 
ity got close to ~ —22 mag. The rising time is estimated 
as about 70 days (hereafter, days are in the rest frame). 
After reaching the peak luminosity, the LC declined slowly 
(~ 0.02 mag day~^) for ~ 120 days and then the LC stayed 
almost constant for ~ 20 days until SN 2006gy hid behind 
the Sun. Near-infrared (NIR) LCs in these early epochs are 
consistent with the blackbody temperature obtained from 
the optical spectra and no significant excess was detected 
l|Miller et al]|2010l). No X -ray was detected when the LC is 
rising (|Ofek et al.ll2007l ) but weak X-rays may have been de- 
tected during the declining phase jSmith et al 1 120071'). No ra- 
dio emissio n is d etected at any observed epochs (lOfek et al.l 
20071: IChandra. Chakraborti. fc Ravll2007l : lArgo et al.ll2007l : 



Bietenholz fc Barte]||2007l ■ l2008al lb^ 



Optica l spectra of the ear l y epochs are also take n inten- 
sively fe.g.. lSmith et al.ll2010l : lAgnoletto et al.ll2009l ). Spec- 
tra taken before the LC peak are charac terized by Lorentzian 
Ha emission lines l|Smith et al.ll2010l ). The origin of the 
Lorentzian profil e is related t o the existence of optically 
thick CSM fe.g.. lChugaill200ll : lDessart et al.ll2009l ). Except 
for the narrow H emission lines, the spectra are feature- 
less and characteri z ed by blackbody with depletion in blue 
jSmith et al.l I2OI0I : lAgnoletto et all boogl ). The reason for 
the lack of features may be partly because the spectra be- 
fore the LC peak were taken only with low resolutions. After 
the LC peak, overall Hq line profiles can be fitted by two 
Gaussian components with the FWHM of 1,800 km and 



5,200 km s^^ and they are presumed to come from the in- 
teracting region between the ejecta and dense CSM (e.g.. 
Smith et al.ji2010i ). There also exists a broad absorption in 
the blue part of the Ha profile which is suggested to origi- 
nate from the ejecta inside. In addition, the spectra of SN 
2006gy show narrow P-Cygni profiles from several elements 
(e.g., H, Fe) with the outfiowing velocity ~ 100 km s~^. As 
the velocity is too slow to attribute it to the ejecta inside, 
those narrow lines are presumed to originate from the un- 
shocked CSM. The strengths of these narrow lines decline 
with time and they are barely se en in the spectrum taken 
at ~ 140 days since the LC peak (|Smith et al.ll2010l ). 

About 100 days later, SN 2006gy came out of 



the Sun and was observed again jAgnoletto et al.l l2009l : 
iKawabata et al.ll2009l : [Miller et al.ll2010l ). The optical lumi- 
nosity of SN 2006gy declined dramatically (about 2 mag 
in the R band) which was almost constant for ~ 20 days 
before the SN went behind the Sun. The luminosity de- 
clined very slowly (~ 0.002 mag day^^) since it appeared 
from the Sun for more than 400 days until the last reported 
observation on Novemb er 22, 2008 (UT) (Kaw abata et al.l 
I2OO9I : iMiller et al1l201Gl '). The dechne rate is much slower 
than that of ^®Co decay (0.01 mag day~^) and the main 
source of the luminosity cannot be the ^^Co decay. Be - 
cause of th e high NIR luminosities, ISmith et al.l l|2008l ): 
I Miller et al.l l|2010t ) suggest that the late time luminos- 
ity is due to light echoes. The optical spectra of those 
epochs are domin ated by intermediate w idth emission lines 
(~ 2,000 km s~\ IKawabata et all [20091 ') . H emission lines 
were weaker than those observed in previous epochs and 
suggest that the interaction is weak in tho se epochs and it 
is no longer a main source o f the radiation ([Agnoletto et al.l 
I2OO9I : IKawabata et al.[|2009l ). Weakness of Fe lines in those 
epochs seems inconsis tent with the large ^®Ni production 
(IKawabata et al.ll2009l ). 



3 NUMERICAL METHOD 

We briefly summarise the basics of STELLA code, which we 
use for our LC modeling. We then describe in detail a key 
numerical method of STELLA that treats smearing due to 
multidimensional effect. 



3.1 STELLA 

STELLA is a one-dimension al multigroup radiati o n hy- 
drodynamics code (e.g., iBlinnikov fc Bartunovl 19931 : 
iBlinnikov et al.ll 19981 . [20061 : IBlinnikov fc Tolstovll201ll ') and 
calculates the spectral energy distributions (SEDs) at each 
time step. Multicolor LCs can be obtained by convolving fil- 
ter functions to the SEDs. All the calculations are performed 
by adopting 100 frequency bins from 1 A to 5 x 10'*A in 
log scale. STELLA implicitly treats time-dependent equa- 
tions of the angular moments of intensity averaged over 
a frequency bin. Local thermodynamic equilibrium is 
assumed to determine the ionization levels of materials. 
STELLA has been intensively used for m odeling the SN LCs 



powered b y the sh ock interaction ( e. g., Chugai et al. 2004 : 
Wooslev. B linnikov fc Hegeij [20071 : IBlinnikov fc Sorokinal 
2010l:lMoriva et aLlbOllI) as well as other types of SNe (e.g. , 



Baklanov. Bhnnikov. fc PavlvukI [ioosi : [FolatelU et al] [20061 : 
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iTsvetkov etHI l2009l : iTominaga eFall I2OIII ). Comparisons 
of STELLA with other numeric al cod e s are provided i n, e.g. , 
iBhnnikov" et all |l99a. ,200o|. I2OO3I ): IWooslev et aD ()2007l l 
and analytical models are also compared to the numerical 
results of STELLA (e.g.. iRabinak fc WaxmanlbOHl ). 



3.2 Smearing 

As is shown i n th e prev i ous works (e.g.. 
Chevalier fc FranssonI 1 19941 : IChugail I2001I : IChugai et al.l 



2004 ) and will be shown in the following sections, the inter- 



action of SN ejecta and dense CSM results in a dense cool 
shell between SN ejecta and dense CSM. This is because 
of the radiative cooling of the shocked region. The shocked 
region becomes very dense and cools down efficiently by 
radiation. The cooling prevents the pressure from growing 
sufficiently enough to sustain the shell. Hence, the shell 
becomes thinner and denser and the cooling becomes 
more efficient. Thus, this cooling process is catastrophic. 
However, in reality, such a shell is unstable because of 
several instabilities like the Ray leigh- Taylor instability 
which require multidimens ional calculations to treat (e.g., 
IChevaher fc BlondinI 1 19951 ). The multidimensional effect 
smears the shell and less kinetic energy is converted to 
radiation energy. In other words, the cooling by radiation 
is less efficient in three dimensions than in one dimension. 
In STELLA code, we take such multidimensional effects into 
account by introducing a smearing term in the equation of 
motion so that the conversion efficie ncy from kinetic energ y 
to radiation energy can be reduced jBlinnikov et al.lll998h . 
This term is similar to artificial viscosity, although the 
smearing term has c ompletely opposi t e effec t. 

As is shown in iBlinnikov et al.l (| 19981 ). the smearing 
term is defined such that the total energy is manifestly 
conserved. Only the neighbouring zones are affected by the 
smearing term. The overall normalization factor i?cut(T) of 
the smearing is expressed as 



7?cut(T) = BJ{t). 



(1) 



See IBlinnikov et al.l (|l998l ) for the definitions of 7?cut(T) and 
the smearing term, /(r) is introduced so that the artifi- 
cial smearing is reduced at optically thick regions where the 
effect of cooling is less efficient and the multidimensional 
instabilities due to the cooling grow less. In STELLA code, 
/(r) is an empirically obtained monotonically decreasing 
function which satisfies /(r ^ 0) = 1. Bq determines the 
overall strength of the smearing effect. Ideally, Bq should 
be calibrated by comparing results obtained by our one- 
dimensional calculations to those of multidimensional calcu- 
lations in which the effect of multidimensional instabilities 
in the shell are taken into account. However, such a multi- 
dimensional radiation hydrodynamical code with which we 
can compare our results is not available yet. We use Bq — 1 
as our standard value. We also show the effect of Bq on 
model LCs (Section 15.2.41 and Appendix \K\ . The results of 
LC calculations strongly depend on Bq , as the smearing term 
directly affects the conversion efficiency from kinetic energy 
to radiation. We need a multidimensional radiation hydro- 
dynamics code for the calibration of the parameter Bq. 



4 INITIAL CONDITIONS 

In this section, we show how the initial conditions for our 
LC calculations are constructed. Two components exist in 
the initial conditions: SN ejecta inside and CSM outside. We 
assume that there is a gap between the progenitor and the 
dense CSM. Then, it takes some time for the SN ejecta to 
reach the dense CSM and start to collide. We assume that 
the SN ejecta freely expands in the gap before the collision. 
We numerically follow the LCs after the collision. The ini- 
tial conditions of the two components, SN ejecta and dense 
CSM, are constructed by the way explained in this section. 
Then, these initial conditions will be confirmed by our nu- 
merical LC calculations as we show in the later sections. 
The initial density structures of two representative mod- 
els are shown in Figure [T] as examples. Both SN ejecta and 
CSM are assumed to have solar metallicity and no *^Ni is 
included in our calculations unless it is otherwise mentioned. 
The summary of the models is given in Table [T] 



4.1 Supernova Ejecta 

SN ejecta before the collision is assumed to be freely ex- 
panding with a homologous velocity profile. The analytic 
approximatio n for the density st ructure of SN ejecta pro- 
vided by, e.g.. lChevalier fc Sokeij (jl989 ) is adopted: 



p{r,t) 



(2) 



where Mcj is the SN ejecta mass, t is time since the explosion, 
vt = Cv {Eb]/ Mf,{f'^'^ , and -Ecj is the kinetic energy of the 
SN ejecta. The constants C,p and are constrained by the 
condition that the sum of density and kinetic energy should 



be Mej and E^f 



Cp 



2(5-5) (n-5) 



1 1/2 



. (3- 
1 (n 



5) (n - 3) 
- 3) (3 - 6) 



4n 



(3) 
(4) 



In this paper, we adopt S — 1 and n = 7, which is used by 
IChevalier fc IrwirJ (|201ll ) . We have also tried n = 6 and n = 
8 for some models but results had little difference compared 
to those of 71 = 7. The maximum velocity of the SN ejecta 
before the interaction is chosen to be high enough, so that 
most of the assumed Eej is contained in the SN ejecta. It is 
around 20, 000 - 50, 000 km s"^ 

A/cj is difficult to be constrained only by the observa- 
tions of the LC of SN 2006gy because the LC is mainly af- 
fected by CSM as is shown in the following sections. In most 
of the models, we adopt Mcj — 20 Mq because the progen- 
itors of Type Iln SNe are presumed to be originated from 
relatively massive stars. Effect of Moj on LCs is discussed in 
Sections 15.2.31 and 16.11 



4.2 Circumstellar Medium 

Dense CSM in the calculations is assumed to exist from Ri 
to Ro from the center and have a density structure p (x . 
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Table 1. List of LC models 



Name 








w 


Ro 


yiRo 


xRo 


Ri 










km s^^ 


Mq 


10^1 erg 




10^5 cm 10^5 cm 


cm 


10-"^^ em 


Mq 


Mq yr-i 




Al 


5,200 


20 


50 


5 


11 


4.9 


1.8 


1.0 


22 


0.70 


1 


A2 


5,200 


20 


50 


5 


11 


4.9 


1.8 


1.5 


10 


0.33 


1 


Bl 


5,200 


20 


10 


2 


8.8 


3.3 


0.090 


0.090 


0.83 


0.030 


1 


B2 


5,200 


20 


30 


2 


8.8 


8.3 


2.0 


0.090 


27 


0.98 


1 


CI 


5,200 


20 


50 





4.9 


4.8 


1.8 


1.0 


14 


1.1 


1 


Dl 


10 000 


20 


10 


5 


21 


11 


4.6 


4.5 


22 


0.42 


1 


D2 


10,000 


20 


10 


5 


21 


11 


4.6 


5.0 


18 


0.36 


1 


El 


10,000 


20 


30 


2 


17 


6.3 


0.32 


0.10 


3.2 


0.060 


1 


E2 


10,000 


2 


10 


2 


17 


6.3 


0.32 


0.10 


3.2 


0.060 


1 


Fl 


10,000 


20 


10 





11 


11 


4.6 


5.0 


15 


0.45 


1 


D3 


10,000 


10 


10 


5 


21 


11 


4.6 


5.0 


18 


0.36 


1 


D4 


10,000 


30 


10 


5 


21 


11 


4.6 


5.0 


18 


0.36 


1 


D5 


10,000 


10 


5 


5 


21 


11 


4.6 


5.0 


18 


0.36 


1 


D6 


10,000 


20 


10 


5 


21 


11 


4.6 


5.0 


18 


0.36 


0.33 


D7 


10,000 


20 


10 


5 


21 


11 


4.6 


5.0 


18 


0.36 


3 


Average mass-loss rate. 


(m) = 


AfcSM 


/{Ro- 


Ri) where Vn, 


= 100 km 


s -"^ is the assumed CSM velocity. 





w — 2 density structures and most of the CSM is likely from 
non-steady mass loss (e.g., [Pw arkadas fc Gruszko 201 li ). A 
steep CS M density gradient wi t h w = 5 is suggested for SN 
2006gy in iMoriva fc Tominagal (|2012l '). We also show w = 
models which are difficult to be excluded only by the LC 
modeling. 

It turns out in the later sections that it is difficult to 
estimate Vs from observational values self-consistently. This 
is partly because Vs is not an independent parameter and, 
in principle, can be derived for a given CSM structure if we 
specify Ecj and Moj. However, Va is also strongly affected 
by the conversion efficiency from the kinetic energy to ra- 
diation through the interaction and it is unknown at first. 
Thus, it is difficult to estimate Vs from the first principles. 
Hence, we set Vs as a free parameter in this paper. At first, 
we try to estimate it from the observations. As the black- 
body radius of SN 2006gy expands linearly with the velocity 
5,200 km s"'^, one may estimate that Vs — 5,200 km s~^. 
However, the required Ecj for the Vs = 5, 200 km s^^ mod- 
els to explain the peak luminosity of SN 2006gy is found 
to be very high and it becomes inconsistent with the rela- 
tively low Vs ■ In other words, the SLSN models obtained by 
setting Vs — 5,200 km s~^ are not self-consistent. Thus, we 
also try models with higher Vs, namely, Vs = 10, 000 km s~^. 
The Vs = 10,000 km s~^ models are found to work well 
self-consistently as is shown in the following sections. In ad- 
dition, the linear evolution of the blackbody radius with 
5,200 km s~^ is found to be able to be explained by the 
Vs = 10,000 km s"^ models. With td = 70 days, ts = 194 
days, and the given lis, we can derive Ro and xRo from 
Equations ((5)1 and (|6]) for a specified w. In the rest of this 
section, we show the details of the two Vs models. 

4.2.1 Vs = 5,200 km s~^ Model 

loss and they are naturally expected structures for CSM. This model corresp o nds to the SN 2006gy model in 
Note that X-ray observations of Type Iln SNe suggest that iMoriva fc Tominagal (|2012l ). The shock velocity Vs = 
CSM around Type Iln SN progenitors often does not have 5,200 km s~^ is estimated from the observed evolution of 



The outfiowing velocity of the CSM is 100 km s ^. It is es- 
timated from the narrow P-C ygni profile of Hq appeared in 
the spectra of SN 2006gy (e-g.- ISmith et al.ll2010l ). The CSM 
is assumed to be optically thick enough to cause the shock 
breakout within it. We estimate the physical conditions of 
the CSM from the observ ations by using the shock b reak- 
out condition described in IMoriva fc Tominagal (|2012l ). Our 
main purpose of this paper is to see how well the properties 
of the dense CSM predicted by this simple shock breakout 
model in CSM can explain the overall LC features of SN 
2006gy. We use following three values which can be esti- 
mated from the observations to derive the CSM properties: 
the photon diffusion time td in CSM, the propagation time 
ts of the forward shock through CSM, and the forward shock 
velocity Vs . td corresponds to the rising time of the LC and 
ts corresponds to the time when the narrow P-Cygni Hq 
profiles from C SM disappears. We adopt td = 70 days and 
ts = 194 days (|Moriva fc Tominagal [20l3 ) . Vs can be esti- 
mated from the spectral evolution. 

With td, ts, and Vs, we can estimate the outer radius Ro 
of the CSM and the radius xRo where the shock breakout 
occurs {Ri/Ro < a; < 1) for a given w based on the shock 
breakout model. The shock breakout c ondition predicts the 
follow ing relations for the three values (|Moriva fc Tominagal 
120121 ): 



td 



{W = l), 



(5) 



(6) 



We try three values for w. w — 0, 2, 5. The models 
with w — 2 corresponds to the case of the steady mass 
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log radius (cm) 

Figure 1. The density structures of the models D2 (ui = 5) and 
Fl (ui = 0) before the coUision. These density structures are set 
as the initial conditions of the numerical LC calculations. The 
structure inside the density jump corresponds to the SN ejecta 
which is assumed to be freely expanding before the collision (Sec- 
tion [4TTJ. The dense CSM exists above the density jump. The 
initial condition of the dense CSM is obtained based on the shock 
breakout model within the CSM (Section l4.2l l. 



the blackbody radius of SN 2006gy (ISmith et aDl2010l ). For 
10 = 5, we get Ro — 1.1 x 10^® cm and xRo = 1.8 x 10^® cm 
from Equations ((5]) and ([6]). The inner radius Ri of CSM can- 
not be constrained by the above observables and it is a free 
parameter. Two Ri are tried: 10^^ cm (Al) and 1.5x10^^ cm 
(A2). With u) = 2 (steady mass loss), we get Ro = 8.8 x 10^^ 
cm and xRo — 9.0 x 10^^ cm (Bl). We set Ri = xRo in the 
w = 2 models. The results do not depend so much on Ri 
in this case because most of the mass in the w = 2 CSM 
is distributed in the outer part of the CSM. We also calcu- 
late LCs from the model (B2) in which the CSM mass is 
artificially increased 30 times of the model Bl. 

In w = models, ts be comes similar to td {ts — 
U) (|Moriva fc Tominagal [2013 ). Thus, the entire CSM is 
shocked with td and no unshocked wind remains after the 
LC peak. This is against the observations of SN 2006gy be- 
cause narrow P-Cygni profiles are observed after the LC 
peak. However, this is only true when we only think a single 
w = CSM component. If there is another CSM compo- 
nent outside the main CSM which is not dense enough to 
affect the LC but the spectra, SN 2006gy-like SNe can ap- 
pear. Thus, w = Q models are difficult to be excluded only 
by the LC. From the rising time of the LC, we can pre- 
sume yiRo — xRo Ro — xRo ~ Vatd — 3.1 x 10^^ cm. 
yiRo is the radius where the optical depth from the surface 
of the CSM becomes 1 an d yiRo — Ro when w = Q (see 
iMoriva k, Tominagal l20lj for the details). We set the last 
scattering surface of the to = 5 model {yiRo = 4.9 x 10^^ cm) 
as Ro so that we can compare the results with those of the 
w = 5 models. Thus, xRo = 1.8 x 10^^ cm is also the same 
as the w = 5 model and we adopt Ri — W^"' cm (CI). 

4.2.2 Vs = 10,000 km s"^ Model 

Vb = 10, 000 km s~^ models are constructed by following the 
same way as the Vs = 5, 200 km s~^ models. With td = 70 



days, ts = 194 days, and Us = 10,000 km s~^, Ro = 2.1 x 
10^® cm and xRo = 4.6 x 10^^ cm are obtained for w — 5. 
We try two Ri: 4.5 x lO^'^ cm (Dl) and 5 x 10^^ cm (D2 
- D7). Although Ri = 5 X 10^^ cm is slightly larger than 
the shock breakout radius xRo = 4.6 x 10^^ cm, it turns 
out that the model gets closer to the SN 2006gy LC. Given 
the approximated way of our estimations, the difference is 
within an acceptable range. The steady mass-loss models 
{w = 2) gives Ro = 1.7 X 10^^ cm and xRo = 3.2 x 10^'' cm 
(El, E2). For w = 0, Ro ~ xRo ~ Vstd = 6 x 10^^ cm (c.f. 
yiRo — Ro). As yiRo ~ 1.1 x 10^^ cm in ui = 5 model, we 
adopted Ro — 1.1 x lO^'' cm and xRo — Ri — 5 x 10^^ cm 
(Fl). 



5 LIGHT CURVE MODELS 

Starting from the initial conditions obtained in Section |4j 
we perform the numerical LC calculations with STELLA. We 
show that the Vs — 5, 200 km s~^ models cannot explain the 
huge luminosity of SN 2006gy self-consistently. LCs obtained 
from Vs = 10, 000 km s~^ models are broadly consistent with 
the observational properties of SN 2006gy. 

5.1 Vs = 5, 200 km s"^ Model 

Figure [5] shows the LCs from the Vs = 5, 200 km s~^ models 
with the observed R band LC of ISmith et"all |20o3). Eej is 
chosen so that the peak luminosities of the model LCs can 
be as luminous as that of SN 2006gy. However, i?ej should be 
very high (~ 5 x lO^'^ erg) for the Vs — 5, 200 km s~^ models 
to be as luminous as SN 2006gy and assuming the relatively 
low lis = 5, 200 km s~^ is not consistent with the high kinetic 
energy. This inconsistency can also be seen from the rising 
times of the models. Although the models are constructed 
so that the rising times of the LCs become td = 70 days, the 
rising times of the numerical results are much shorter than 
70 days. If we set smaller _Ecj, the rising times can be the 
same as that of SN 2006gy but then the luminosities become 
much smaller than that of SN 2006gy. Shortly, the models 
derived by assuming Vs = 5, 200 km s~^ are not compatible 
with the large luminosity of SN 2006gy. Given these results, 
we adopt models with a higher Vs, Vs = 10, 000 km s~^, and 
they are found to be able to explain the LC of SN 2006gy 
self-consistently (Section 15. 2|l . 

One important question of the interaction model is 
whether the CSM from the steady mass loss (w = 2) can 
explain the properties of SLSNe and we look into the w — 2 
models more carefully. If it can, a mechanism to achieve such 
huge steady mass loss may exist. If not, it is indicated that 
explosive non-steady mass loss takes place in their progeni- 
tors and there should exist some mechanisms to cause such 
mass loss just before their explosions. 

The w = 2 model Bl reaches only ~ —19.5 mag in the 
R band at the LC peak. This is because Mcsm = O.83M0 
in the Bl model is much smaller than Mcsm of the models 
with the other w. The fraction of kinetic energy converted 
to radiation in the model Bl is much smaller than those in 
the models Al, A2, and CI, as the amount of energy con- 
verted from kinetic energy to radiation strongly depends on 
the relative mass of CSM and SN ejecta (see Section [6.ip . 
Thus, more kinetic energy is required for the Bl model to be 
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Figure 2. Absolute R band LCs of the Vs = 5,200 km 8"^ 
models. These models are not self-consistent. The origin of time 
axis is set to when our numerical calculations start, i.e., when the 
SN ejecta and CSM start to collide, except for B2. The time of 
the model B2 is shifted —20 days. 
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D1 (w=5, Ri=4.5x10'f^cm) [-10 days 
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E2 (w=2, Mei=2 Msun, Eei=10 erg 
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Figure 3. Absolute R band LCs of the Vs = 10,000 km s"! 
models. These models can self-consistently explain the rising time 
and the peak luminosity of SN 2006gy. The origin of the time axis 
is -1-10 days (Dl), +5 days (D2), days (El, E2, and Fl) since 
the collision. 



as luminous as SN 2006gy. However, the rising time of the 
Bl model is already much less than that of SN 2006gy and 
it becomes shorter if we increase kinetic energy. Thus, the 
Vs = 5, 200 km s~^ model with the steady wind (w — 2, Bl) 
is hard to be compatible with SN 2006gy. For demonstration, 
we also calculated a model (B2) in which Mcsm is increased 
30 times more than that of the model Bl. Then, the amount 
of energy converted increases because of the high efficiency 
for the energy conversion. In addition, the photospheric ra- 
dius is increased due to the increased density. As a result, 
the luminosity of the model becomes as large as that of SN 
2006gy. 



As lis = 5, 200 km s ^ models are not able to explain SN 
2006gy self-consistently, we investigate models with higher 



= 10, 000 km s" 



— 5, 200 km s is estimated 



from the evolution of the blackbody radius but it is shown 
that the evolution of the blackbody radius in the Vs — 
10, 000 km s"'^ models is consistent with that of SN 2006gy. 



5.2.1 Light Curve 

The R band LCs from the Vs = 10, 000 km s~^ models are 
shown in Figure [S] Multicolor LCs of the models D2 {w — 5) 
and Fl {w = 0) are shown in Figure 3] and the bolometric 
LCs of the two models are shown in Figure [S] The color 
evolution of the models D2 and Fl are shown in Figure [6] 

The rising parts and the peak luminosities of the LCs 
of the w — 0,5 models are consistent with SN 2006gy. Thus, 



the Vs = 10, 000 km i 



models are self-consistent with the 



assumed iSoj and A/cj. The Vs = 10, 000 km s~^ models only 
require iSoj ~ 10^^ erg to achieve the peak luminosity of SN 
2006gy, instead of iSej — 5 x 10^^ erg required for the Vg = 
5, 200 km s"^ models. This is because the blackbody radius 
in the CSM can be larger in the Vg = 10, 000 km s~^ models 
and less energy is required to achieve the same luminosity. 
The steady mass-loss models {w = 2) are, however, still not 
consistent with SN 2006gy. 

The model LCs with w — 0,5 after the peak start to 
deviate from the observed LCs, although the deviations stay 
less than one magnitude before the plateau in the observed 
LC at around 200 days. Our model R band LCs take some 
time after the LCs have reached the peak until the LCs 
start to decline, contrary to the observed R band LC of 
SN 2006gy. This is because there remains unshocked opti- 
cally thick CSM even after the LC peak in our numerical 
models and the p hotosphere remains t h ere fo r a while. The 
analytic model of lMoriva fc T ominaga' (|2012l ) which we use 
for the estimate of the initial condition assumes a constant 
Vs. However, Vs actually reduces as the interaction goes on 
and the optically thick part of CSM does not shocked away 
entirely at the time td when the optically thick CSM is as- 
su med to be swept up by the forward shock in the model 
of iMoriva fc Tominagal l|2012l ). This effect is more signifi- 
cant in w = models because to = models suffer more on 
the deceleration than w = 5 models. A severe failure of our 
models is that all of them fail to reproduce the plateau in 
the LC of SN 2006gy at around 200 days. We discuss this 
separately in Section [6.21 

Looking at the multicolor LCs (Figure \^ , the w = 
model (Fl) is closer to the observed LC, especially the B 
band LC of the rising epochs. This is presumed to be be- 
cause the initial density jump between SN ejecta and CSM 
is smaller in the w = model (Figure [T} and the tempera- 
ture becomes lower in the w = model. However, the LCs 
in the U and B bands can be affected by many weak ab- 
sorption lines of Fe group elements which are not taken into 
account in our opacity. Those weak absorptions may reduce 
the luminosity of the U and B band LCs and we cannot 
judge which model is better based just on the blue part of 
the LCs. In addition, the difference is ~ 0.5 mag and they 
are not significant. The color evolution of the two models 
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(Figure [6]) roughly follows the observed evolution, although 
there exist some deviations especially in R — I and V ~ R. 

Although it is possible to continue LC modeling to get 
much better fits to the SN 2006gy LC, it does not neces- 
sarily lead us to the better understanding of the proper- 
ties of the SN ejecta and the dense CSM involved in the 
progenitor system of SN 2006gy. This is mainly because of 
the simplified physics adopted in STELLA. Especially, STELLA 
is a one-dimensional code and multidimensional effects are 
approximately incorporated by adopting the smearing pa- 
rameter. As the uncertainties involved in the parameter are 
large (see Section [5]2]4]), making a perfect fit to the observed 
LC does not necessarily provide us with the best parame- 
ter. In addition, the differences in the LCs in the declining 
phases are less than one magnitude (or a factor ~ 2) and 
the differences in the rising phases are much less. Thus, the 
properties of the SN ejecta and the dense CSM in the D2 
and Fl models are presumed not to be so different from the 
'actual' values. Thus, we conclude that the CSM parameters 
predicted by the shock breakout model can explain the over- 
all properties of SN 2006gy. We also note that a systematical 
study of the effect of the CSM properties on the LCs pow- 
ered by the inter action between SN ejecta and dense CSM 
is summarized in lMoriva et al.l l|201l[ ). The durations of the 
LCs of the models D2 and Fl are a bit longer than that of 
SN 2006gy. To reduce the durations of the LCs by keeping 
the peak luminosities of them, we can, for example, change 
the radu of the CSM. 

Finally, we look into the steady mass loss models to 
see whether they actually fail to reproduce the LC of SN 
2006gy. The steady mass-loss model (w = 2, El) is, again. 



The rising time is already too short and reaching the peak 
luminosity of SN 2006gy by increasing i?ej does not work as 
is discussed in Section [5.11 This is because of the too small 
McsM and can be improved if AfcsM is increased (see the 
models Bl and B2). Another possible way to make w = 2 
models work is to increase the conversion efficiency from iJcj 
to radiation energy so that i?ej can be reduced (see Section 
16.11 for the discussion of the conversion efficiency). In the 
model E2, Mcj is set to be comparable to Mqsm so that 
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Figure 5. Bolometric LCs of the models D2 (ui = 5) and Fl 
(to = 0). The origin of the time axis in the D2 (Fl) model is 5 
(0) days since the collision. 



the conversion efficiency becomes higher (Section 16. Ijl . The 
similar peak luminosity to the model El is reached with less 
i7ej (lO^'^ erg) in the model E2. However, as the diffusion 
time of the CSM is not affected so much by this, the LCs 
become similar to each other and increasing the efficiency 
does not revive the w = 2 models. To summarize, the dense 
CSM from the steady mass loss is still difficult to explain 
the LC of SN 2006gy with the shock breakout model. 



5.2.2 Dynamical Evolution 

Figure [7] shows the dynamical structures of the model D2 
{w = 5) at around the LC peak. The left panel shows the 
structure in the physical coordinate (radius) whereas the 
right panel shows the structure in the mass coordinate. The 
cool dense shell is created between SN ejecta and CSM in 
which about 10 Mq of shocked CSM and 3 Mq of shocked 
SN ejecta are contained. 

Figure |8] shows the evolution of the color temperature 
(Tcoi) and the effective temperature (Tcff) of the model 
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Figure 7. Physical structures of the model D2 in radius (left) and mass coordinate (right) at around the LC peak (74 days since the 
collision). Black lines show the density structure (left j/-axis). Blue dotted lines are the velocity scaled by 10* cm s~^ (right j/-axis), 
purple lines are the logarithm of the absolute value of luminosity scaled by lO**" erg s~^ (right y-axis), green lines are the logarithm of 
the temperature in Kelvin (right y-axis), and red lines are Rosseland optical depth measured from the outside (right y-axis). 



D2. The color temperature is derived by fitting the spec- 
tra obtained by the numerical calculations with the black- 
body spectral distribution whereas the effective temperature 
is obtained by using the bolometric luminosity (Lboi) and 
the radius (i?^^^2/3) of the photosphere which is defined 
as the radius where the Rosseland mean optical depth th 
from the surface becomes 2/3 in STELLA and is expressed 

as Tcff — (i'boi/47rcr_R^^^2/3)^^''- Here, cr is the Stephan- 
Boltzmann constant. As radiation mainly comes from the 
shell and the Thomson scattering is the dominant opac- 
ity source in the CSM above the shell, Tcoi roughly traces 
the temperature of the shell. The photosphere (-Rtj,=2/3) is 
much above the shell and Tefi becomes very low because of 
the large -Rtjj=2/3 (see also Figurel^)). At the time when Tcoi 



starts to increase for the second time (from ~ 180 days), the 
photosphere is in the SN ejecta whose density structure and 
composition are expressed in the approximated way and the 
results around these epochs and later should not be taken 
seriously. 

Figure [9] shows the evolution of the blackbody radius 
-Rbb ~ \J I/boi I'^T^^TtoX photosphere -Rtj,=2/3- The 

constant velocity line with 5, 200 km s~^ is the ev o lution 
of the blackbody radius obtained by ISmith et al.l (|2010l ) 
along which 7?bb follows until around 125 days. The evo- 
lution of /?BB in the models D2 and Fl is consistent with 
5, 200 km s^^ although the radius is a bit sma ller than the 
observed values. The bolometric correction of ISmith et al.l 
((2010.) is based on Tbb and the correction may add extra 
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Figure 8. Evolution of thie color temperature (Tcoi) and the 
effective temperature (Tgff) in the model D2. The origin of the 
time axis is 5 days since the collision. 



luminosities because it ignores the effect of line depletion. 
Higher Lboi results in higher Rbb for a given Tco i and this 
can be the reason why 7?bb of lSmith et al.l l|201Cf ) is higher 
than ours. 

-Rbb obtained by our calculations tends to be smaller 
than the shell radius where the radiation is coming from. 
For example, i?BB of the model D2 shown in Figure [5] 
stays lower than the radius at which the interaction starts 
{Ri = 5 X 10^^ cm). The reason is presumed to be similar 
to that of the discrepancy in -Rbb obtained from observa- 
tions and numerical calculations. Tcoi is obtained from the 
spectral fitting but actual spectra suffer from the line de- 
pletion especially in blue. Since Ted is reduced from the 
temperature at the photon production site, Lboi is the value 
affected by such depletion and is less than the value ex- 
pected from the blackbody with Tcoi. Thus, with the smaller 
^boi, Rbb ~ \/ -i/boi/47raT*jj[ becomes smaller than the ac- 
tual emitting region. 

As most of hydrogen in CSM remains to be ionized, 
^Tj,=2/3 continues to be at the radius where the Rosseland 
mean opacity from the surface of the CSM is 2/3 and re- 
mains to be constant until the shock wave comes close to 
the radius. Then, -Rth=2/3 evolves roughly foUowing the for- 
ward shock. 



The properties of SN ejecta (-Mej and -Eej) determine many 
aspects of SNe powered by the shock interaction (e.g., lumi- 
nosities) because -Eej determines the available energy and 
-Moj affects the efficiency to convert the available kinetic 
energy to the radiation energy. We discuss the effect of 
-Mej and -Boj in Section 16.11 including the results of the 
lis = 5, 200 km s^^ models and here we just show the re- 
sults of LC calculations with different A/cj and -Ecj (Figure 
nop . LCs are similar to each others and we can see that it 
is difficult to constrain Mej and .Ecj only by the LC. This 
can also be seen by comparing the models El and E2 in Fig- 
ure [S] The two models have different -Mej and -Bej with the 
same CSM but the resulting LCs are similar (see discussion 
in Section [OT) . 
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Figure 9. Evolution of the blackbody radius (-Rbb, solid line) 
and the photospheric radius (.R.r_2/3, dashed line) of the mod- 
els D2 {w = 5) and Fl {w = 0). The origin of the time axis 
is 5 days since the collision. The monotonically increasing linear 
line at th e middle is the evo lution of the blackbody radius ob- 
tained bv lSmith et al. a constant velocity evolution with 
5,200 km s~^. The observational blackbody radiu s follows the 
line until around 125 days and starts to decline fsee lSmith et all 
{2OIQ) for details). 
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Figure 10. R band LCs from the same CSM but different M^j 
and -Eej. The models D2, D3, D4 have the same E^j but different 
-Mej, i.e., 10 Mq (D3), 20 Mq (D2), and 30 M© (D3). The model 
D5 has less Egj (5 X 10^^ erg) than other models (10^^ erg) and 
-Mgj = 10 Mq. The calculation of the LC of the model D3 stopped 
at around 175 days since the collision and we show the LC until 
around 175 days. 



In the best LC m odel of the pulsational pair-in st abilit y 
model presented by IWooslev. Blinnikov. fc Hegeij |20o3), 
the ejecta with 5.1 Mq and 2.9 x 10^^ erg collides the CSM 
with 24.5 Mq. Our canonical models (D2 and Fl) have 
much higher -Eej (lO^'^ erg). One of the reasons is presumed 
to be the smaller photospheric radius {yiRo = 1.1 x lO^'^ 
cm) in our models. The dense CSM in the pulsational pair- 
instability model extends to about 3 x 10^^ cm and the pho- 
tosphere can be larger than our models. This effect of the 
locations of the photosphere can also be seen in comparison 
to the Vs — 5, 200 km s~^ models. The photospheric radii of 
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them are only y\Ro — 5 x 10^^ cm at most and the required 
energy to achieve the maximum luminosity of SN 2006gy is 
5 X 10^^ erg which is even larger than 10^^ erg required for 
the Vs = 10,000 kms"^ models {yiRo = 1.1 x 10^^ cm). 
This shows the difficulties to constrain E^.^ only by the LC. 

In addition, the efficiency to convert the kinetic energy 
to radiation is mainly determined by the relative mass of the 
ejecta and the collided CSM. ft does not depend strongly 
on the ejecta mass if the CSM mass is much larger than 
the ejecta mass. To get high conversion efficiencies of the 
kinetic energy to the radiation energy, Afoj is better to be 
comparable or less than Mcsm and we can at least get some 
constraint on Mcj from the LC based on the view point of 
the conversion efficiency (Section 16. 3[) . 
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5.2.4 Smearing 

The dense shell which appears between SN ejecta and CSM 
is unstable in multidimension as is discussed in Section [3.21 
As a result of the instabilities, less kinetic energy is expected 
to be converted to radiation because there would be the 
extra multidimensional motions caused by the instabilities. 
To take into account such multidimensional effects in one- 
dimensional code STELLA, we include a smearing term in the 
equation of motion (the parameter Bq, Section [3.2^ . 

Figure [TT] shows the LCs with different values of the 
smearing parameter Bq. With larger Bq, the effect of the 
smearing becomes larger and less kinetic energy is converted 
to radiation. In other words, radiative cooling becomes less 
efficient. The model D2 is calculated with our standard Bq — 
1. The model D6 has Bq = 0.33 and the model D7 has 
Bq = 3. The shape of the LC is different even if we only 
change Bq with a factor 3. We also show the effect of Bq on 
the LCs obtained from the pulsational pair- instability SN 
models of lWooslev. Blinnikov. fc Hegerl (|2007l ') in Appendix 
|X]and show that the effect is not unique to our models. We 
discuss the efficiency in detail in Section [6.11 

The uncertainty in the smearing parameter adds an- 
other difficulty in our estimations of physical parameters of 
the progenitor system. This is one reason why we think that 
making the perfect fitting now does not lead us to the exact 
parameters of the progenitor system. The calibrations for 
the smearing parameter should be done at least. However, 
the rising time and the peak luminosity is not so sensitive 
to the smearing parameter and the parameters of SN ejecta 
and dense CSM we obtain with the current uncertainty are 
presumed to be close to the real ones. 



5.2.5 Effect of 

We have also examined the effect of ^®Ni decay on the LCs. 
Figure [12] shows the results. We include ^^Ni at the center 
of the model D2. If we include ^^Ni, the length of the peak 
is extended due to the extra heat source. The significant 
effect can only be seen when we include ~ 10 Mq of '^^Ni. 
However, the amount of ^^Ni is observatio nally constrained 
to be less than 2.5 Mq (jMiller et al.ll2010l ) and the effect of 
®^Ni is negligible. 



Figure 11. R band LCs with different Bq, i.e., D2 {Bq = 1), 
D6 {Bq = 0.3), and D7 {Bq = 3). The smearing parameter Bq 
changes the conversion efficiency from kinetic energy to radiation. 
With Larger Bq, more kinetic energy remains and less radiation 
energy is emitted. The origin of tlic time axis is 5 days since the 
colHsion. 




100 150 
arbitrary days (rest frarre) 



250 



Figure 12. R band LCs with 56Ni. 0.1 Mq, 1 Mq, and 10 Mq 
of ^^Ni is included at the center of the model D2. Only 10 Mq of 
^^Ni makes significant effect on the LC. 



6 DISCUSSION 

6.1 Conversion Efficiency 

The source of radiation in our LC models is the kinetic en- 
ergy of SN ejecta. The amount of energy converted from 
kinetic energy to radiation can be estimated by the conser- 
vation laws of energy and momentum. If we assume that 
the radiation pressure does not change the dynamics of the 
materials so much, the conservation of momentum requires 



A/colcj '^colcj ~ (Mcolcj -I- MsCSm) Wshcll, 



(7) 



where Mcoioj is the mass of the collided SN ejecta, «coicj is 
the mean velocity of the collided SN ejecta, Mscsm is the 
mass of the shocked CSM, and VshcU is the velocity of the 
dense shell between SN ejecta and CSM. Radiation energy 
-Brad emitted as a result of the interaction can be derived 
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from the conservation of energy 

-Erad ~ 



2 JWcolcj ^colcj 



1 2 
- (A/colcj + MsCSm) Wsholl 



,(8) 



where a is the fraction of kinetic energy converted to radi- 
ation. From Equations (O and ((Sj, 



I Mr. 



If most of the SN ejecta and CSM are shocked, i.e., Mcoioj — 
Mcj and Mscsm — Mcsm, we get the rough estimate for the 
radiation energy emitted 

qMcsm 



' MsCSM 



(9) 



E 



Mcj + McsM 



Er 



(10) 



a is expected to be close to 1 without the smearing parame- 
ter Bq because most of thermal energy gained by the shock 
is eventually emitted as radiation. Since the parameter Bq 
adds additional acceleration to reduce the amount of energy 
converted to thermal energy, a is expected to become lower 
as Bq becomes larger. The rest of energy is mostly in the 
form of kinetic energy. We may also express the effect as 
the reduction of the radiative cooling efficiency because less 
radiation energy is emitted with the smearing term. 

Table [2] is the list of radiation energy which is obtained 
by adding up the bolometric luminosity from the time of 
collision to around 300 days since the collision. The model 
D3 is excluded because we do not have the entire numeri- 
cal LC. We also show the parameter a which is derived by 
using Equation (|10p . The efficiency -EradZ-Eoj to convert SN 
kinetic energy to radiation is plotted in Figure [13] as a func- 
tion of McsM/( - Mei -f- Mqsm) with the results obtained by 



Ivan Marie et all (|2010| ). 

At the high Mcsm/ (Afoj -I- Mcsm) region, our standard 
Bq — 1 results follow the line of a = 0.5. This means 
that the efficiency to convert the kinetic energy to radia- 
tion is reduced by 50%. On the other hand, the results of 



tion IS reaucea py 
Ivan Marie et al.l |2 



2010) follow the a = 1 line and the ef- 
fect of multidimensiona l insta bilities is not significant. Al- 
though [van Marie et all l|2010l ) use a three-dimensional code 
and multidimensional instabilities are included in principle, 
their approximated way to treat the radiation and limited 
spatial resolution may have prevented multidimensional in- 
stabilities from growing. 

As Mcsm/ (Afoj -I- A/csm) becomes lower, the results 
start to deviate from the constant a line. This is because 
Mcsm gets very small and most of the ejecta is not af- 
fected by the interaction, i.e., the assumption Mcoiej = Mcj 
is no longer valid and Equation HlOp should not be used. We 
should use Equation ((S)) instead. As Afcoicj «C Mcj in this 
regime, the efficiencies tend to be higher than the values 
obtained from Equation (|10[) . 

The combinations of Eej and Mej which give a similar 
[Mscsm/ (Mcoiej + Mscsm)] I Mcoiej "ej are expected to result 
in similar LCs and they are degenerated. Thus, it is difficult 
to constrain the exact value for Mcj and Ecj from LCs. This 
is clearly seen in the models El and E2 in Figure [S] Both 
the models have similar LCs. The CSM of the two models 
is exactly the same but and iJej are different. Although 
Mcoiej < Mej in the model El and Mcoiej — Mej in the model 
E2, Mcoioj and iMcoioji'cj are happened to be similar in the 
two models with the similar Mscsm — Mcsm ■ Thus, the two 



Table 2. Conversion Efficiency from Kinetic Energy 
to Radiation Energy 



Name 



lO^i erg 



M„i +McSM 



Al 


13 


0.25 


0.52 


0.48 


A2 


12 


0.24 


0.33 


0.72 


Bl 


0.46 


0.046 


0.040 


1.1 


B2 


8.7 


0.29 


0.57 


0.50 


CI 


14 


0.28 


0.41 


0.68 


Dl 


2.7 


0.27 


0.52 


0.51 


D2 


2.2 


0.22 


0.47 


0.47 


El 


2.6 


0.086 


0.14 


0.63 


E2 


3.0 


0.30 


0.62 


0.48 


Fl 


2.9 


0.29 


0.64 


0.46 


D4 


1.7 


0.17 


0.38 


0.46 


D5 


1.6 


0.32 


0.64 


0.50 


D6 


3.4 


0.34 


0.47 


0.72 


D7 


0.89 


0.089 


0.47 


0.19 



models have similar [Mscsm/ (Mcoiej + Mscsm)] fMcoioj? 
and result in the similar LCs. 



6.2 Origin of the Plateau Phase 

There exists a plateau in the LC of SN 2006gy at around 
200 days. None of our models is succeeded in producing the 
plateau. This is because the remaining CSM at these epochs 
is too thin to affect the LC. Note that the LC observations 
at later epochs reject the possibility to explain this plateau 
by ^®Ni heating (see also Section [5.2.5[) . There are several 
other possible ways to explain the plateau. One possibility is 
the recombination in the SN ejecta. Because we use the sim- 
plified SN ejecta structures, our results of LC calculations 
after the photosphere gets inside of the SN ejecta are be- 
yond the applicability of our simple models. Increasing the 
SN ejecta mass may also help because a plateau phase can be 
longer with larger hydrogen mass, although the conversion 
efficiency from kinetic energy to radiation is also affected at 
the same time fSection 16. ip . By putting more realistic SN 
ejecta with realistic hydrogen-rich envelopes or more massive 
SN ejecta, the recombination wave may stay the envelope for 
a while and may end up with the plateau phase, as i s the 
case in Type IIP SNe (see, e.g.. lKaseri fc Wooslevll2009h . We 
note that the blackbody temperatures of these epochs are 
~ 6000 K and they are consistent with this scenario. Re- 
combination may also occur in the shocked CSM or dense 
cool shell. 

Light echoes from the remaining CSM may also play a 
role. The LC after this plateau phase remains almost con- 
stant for more than 2 00 days, although the luminosity is 
about 10 times smaller (|Miller et al.ll201(3f ). Thus, it is possi- 
ble that there existed another CSM component which caused 
the echoes at around 200 days and shocked away when SN 
2006gy was behind the Sun. 

Smearing may also be relevant to the plateau. In the 
models in Appendix [^ we can see that the plateau can ap- 
pear or dissapear depending on the degree of the smearing. 
Less smearing makes the cool shell denser and photosphere 
can stay there for a longer time, possibly making the plateau 
phase. 
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Figure 13. Efficiency of the c onversion of kinetic ene rgy to radi- 
ation. Most of the results from lvan Marie et ah ll l2010t) . including 
aspherical models, are also shown, a is the measure for the effect 
of the smearing parameter Bq. a = 1 means NO smearing effect 
and the smearing effect increases as a gets small. 
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Figure 14. Observed R band LC and the bolomet- 
ric and R band LCs calculated with th e parameter shown 
in IChatzopoulos. WheelerT fc Vinkol 1 I2OI2I) . The observational 
points are shifted arbitrarily to match the numerical result in 
the figure. 



6.3 Progenitor of SN 2006gy 

As the properties of the progenitor system which can be ob- 
tained from the LC modeling strongly depends on CSM, it 
is difficult to get information on the progenitor from the LC 
of SN 2006gy. The CSM properties and SN ejecta properties 
are degenerated. However, we can get some indications for 
it. As the origin of the luminosity is the kinetic energy of 
SN ejecta, the kinetic energy should be converted to radi- 
ation efficiently. If Mcsm <S Afoj, the conversion efficiency 
Q.Mcsm/ (Afoj + Mcsm) is so small that the kinetic energy 
cannot be converted efficiently enough to explain the LC of 
SN 2006gy. Thus, the mass of the CSM should be close to 
or larger than Mej . 

According to our modeling, Mcsm is required to be 
~ 10 Mq and this means that Moj is expected to be 
~ 10 Mq or less. This indicates that the total mass of 



the system well exceeds 10 Mq and the progenitor of SN 
2006gy should be a very massive star. In addition, the 
progenitor should lose Mcsm within ~ 10 years before the 
explosion. Our models for SN 2006gy have Mcsm — 18 Mq 
and it may be difficult for red supergiants (RSGs) to 
have such mass loss because of the following reason: 
To have CSM with 18 Mq, the zero-age main-sequence 
mass of RSGs should be very large but such massive 
stars suffer more from the radiation driven wind during 
their main-sequence phase because of their large lumi- 
nosities. Thus, losing most of their mass only just before 
their explosions might be difficult. However, extensive 



van Loon et al.l 20051: Vanbeveren, Van Bever. & BelkusI 


20071: ISmith, Hinkle, & Rvde' '2009*: 


'Rover et al.' '20W|; 


Yoon & CanticUo 


2010; Moriya ct al. 


2011; Georgy 2o3) 



and it is still possible that very massive RSGs lose ~ 10 Mq 
just before their explosions due to, e.g . . pulsations (e 



I Yoon fc Caiitielldl2010l:lHeger et al.|[T997l : IlI fc GondllQol. , 
dust (e.g. , van Loon ct all 120051 ) . or g-mode oscillations 
(jQuataert fc Shiodc 201^, see also lArnett fc Mea kin 2011). 

Another possible progenitor of SN 2006gy is a very 
massive star in the luminous blue variable (LBV, see, e.g., 
iHumphrevs fc Davidsonlll994l ) phase w hich is suggested t o 
be the progenitor of SN 2006gy (e.g., ISmith et all 120071 ). 
LBVs experience extensive mass loss an d some of them, 
e.g., 7; Carinae, have ~ 10 Mq CSM (e.g.. [Smith fc Owockil 
I2OO6I ). Their typical wind velocities are also consistent with 
the wind velocities estimated from the narrow P -Cygni Ha 
profiles in SN 2006gy (e.g., ISmith et al.ll2010l ). Although 
LBVs are theoretically considered to be on the way to Wolf- 
Raye t stars and do not explode (e.g.. [Crowtheij [20071 : IVin^ 
[iooi), the progenitor of Type I ln SN 2005gl is f ound to be 
an LBV (|Gal-Yam ct al. 20o3; ICal-Yam fc Leona rd 200^). 
Several other Type Iln SNe are als o suggested to have 
evidences for LBV progenitors (e.g., Kotak fc VinkI 20061 : 
iTrundle et a"l]|2008l : ISmith et al.ll201ll : iKiewe et al.ll2012D 



A shell created due to the interaction between the RSG 
wind and the Wolf-Rayet wind is another possible way to 
have a massive CSM (e.g.. lDwarkadaj|201ll ). Alternatively, 
a shell created by pulsational instability can be followed 
up by the SN ejecta , instead of the ejecta of the nex t 
pulse as suggested bv lWooslev. Blinnikov. fc Hegeij (|2007h . 
Some binary interaction may cau se extensive mass loss (e.g., 
iHachisu. Kato. fc Nomotg I2OO8I ) but binary systems have 
not been considered deeply as a possible progenitor of SLSNe 
yet (see [chevalier 2012; Sokcr 2012). The collision of mas- 
sive stars in a dense stellar cluster can make a massive star 
surrounded by a massi ve CSM and it may also result in 
SN 2 006gy-li ke SLSNe JPortegies Zwart fc van den Heuvell 
I2OO7I . see also lPan. Loeb. fc Kasenll2012l ). 

With the condition that Mej is similar to or less than 
Mcsm, the conversion efficiency of kinetic energy to radia- 
tion (Equation ((5|) is expected to be ~ 50 % at most (Sec- 
tion I6.1|l . As the radiation energy emitted by SN 2006gy 
exceeds 2 x 10^^ erg, the SN ejecta should have more than 
~ 4 X 10^^ erg. Thus, the SN explosion inside should be 
very energetic. As the energy of our models is comparable 
to those of energetic broad-line Type Ic SNe whose progeni- 
tors are suggested to be very massive (e.g., iNomoto et al.l 
the estimated high energy may also indicate that 
the progenitor mass is rather close to those of LBVs. Note, 
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however, t hat the host Rala xy of SN 2006gy is not metal- 
poor (e.g., lOfek et all 120071 ') while broad-line Type Ic SNe 
appe ar more preferentia l ly in low metallicity environments 
(e.R., lArcavi et al.|[2oTol : iModiaz et al.ll201ll : ISanders et"all 
I2OI2I ). In addition, the late time spectra of S N 2006gy are not 
simil ar to those of broad-line Type Ic SNe (iKawabata et al.] 
|2009| ). ahhough the late t ime spectra of Type Ic SLSN 
2010gx show such features (jPastorello et al.ll2010l 'l. 



6.4 Comparison to Semi-Analytic Model 

Recently. IChatzopoulos. Wheeler, fc Vinkd l|2012l ') proposed 
a semi-analytic model of LCs powered by the shock interac- 
tion. The model involves several simplifications but the over- 
all features predicted by the model are shown to match some 
numerical results. They show a LC model of SN 2006gy and 
we perform LC calculations with the same parameters which 
are obtained by them. The parameters are 5 = Q, n = 12, 
w = 0, -Eoj = 4.4 X lO'^^ erg, Moj = 40 Mq , _Ri = 5 x 10^"* cm, 
and Ro — 2.5 x 10^^ cm (corresponds to Mcsm = 5 Mq 
with a constan t CSM density 1.5 x 10~^^ g cm~ ^ ). FoU ow- 



ing the result of lChatzopoulos. Wheeler, fc Vinkol (|2012l ') , we 

put 2 Mq of ^®Ni at the center of the ejecta but the value is 
too small to affect the main part of the LC fSection l5.2.5p . 
Except for the central region, the composition is set as the 
solar metallicity. We use our standard Bq = 1. 

Figure [T?] shows the result of the numerical 
calculati on. Overall, the pa rameters suggested by 
[Chat zo poulos. Wheeler, fc Vinkd 12012) do not result 
in a similar LC to that of SN 2006gy. The peak luminosity 
of the bolometric LC is close to that of SN 2006gy but 
the duration is much shorter that that of SN 2006gy. In 
addition, the R band peak luminosity is much smaller than 
that of SN 2006gy because the photospheric temperature 
is much higher than that of SN 2006gy around the peak 
and the model has much bluer spectra. The duration can 
be longer if we use a smaller Bq but our results shown in 
Figure [11] imply that it is difficult to make the duration 
two times longer than the Bq = 1 LC to match the LC of 
SN 2006gy by just making the Bq small. What is more, 
changing Bq does not improve the color of the LC and the 
R band LC is expected to remain much fainter than the 
observed R band LC. The plateau phase after the drop in 
the LC is due to the recombination in 40 Mq SN ejecta 
inside. The duration of the plateau phase should actually 
be much shorter than that in Figure 1141 as the entire SN 
ejecta is composed of the solar metallicity in the model 
(except for the central ^''Ni) and we use an approximate 
density structure for the SN ejecta. 

There are several possible reasons for the discrepancy. 
The semi-analytic model assumes that the thermal energy 
gained by the forward shock is always released at the center 
of the CSM because the assumption is required to treat the 
transport equations analytically. However, this assumption 
keeps the diffusion time of the photons from the forward 
shock constant and the diffusion time is fixed with the ini- 
tial value. In reality, the forward shock travels outward and 
the diffusion time decreases with time as the remaining un- 
shocked wind decreases. This effect leads to the overestima- 
tion of the duration of the LCs in the semi-analytic model. 
This is presumed to be the main reason why the duration 
of the LC obtained by the numerical calculation is much 



shorter than that obtained by the analytical model. This in- 
dicates that the semi-analytic model should not be applied 
to the system with CSM of a large initial diffusion time in 
which the kinetic energy of shock waves is the main source 
of radiation. 

Another possible reason is that the energy released 
by the reverse shock is overestimated in the semi-analytic 
m odel. In the sem i- analytic model, the self-similar solution 
of lChevalieJ (|l982l ): lNadezhinl l|l985l ) is used as the evolution 
of the hydrodynamical structure. However, in reality, the ef- 
fect of cooling which is not taken into account in the adia- 
batic self-similar solution is so strong in the case of SLSNe 
powered by the shock interaction that a thin cool dense shell 
is created between the SN ejecta and the dense CSM. Thus, 
the reverse shock could not travel as fast as expected from 
the adiabatic self-similar solution and rather stays close to 
the forward shock. 

In summary, many important effects which are essen- 
tial in modeling the LC powered by the shock interaction 
between SN ejecta and a dense CSM with a large photon 
diffusion time lack in the semi-analytic model and it may 
not be appropriate to use it for the modeling of SLSNe pow- 
ered by such a strong interaction. 



6.5 High-Redshift Type Iln Superluminous 
Supernovae 

SLSNe powered by the shock interaction are not only 
brighter but also bluer than other SNe. Thus, SLSNe are 
expected to be a good probe to study the high-redshift 
Universe. Especially, as their progenitors are expected to 
be very massive (Section 16. 3p . SLSNe can provide us with 
the information of very massive stars in the early Universe. 
Some hig h-redshift Type Iln SNe ar e already detected up t o 
2 = 2.36 l|Cooke et al.ll2009l . see also lBarton fc Cookdlioiol ') . 

With upcoming optical deep s urveys with, e.g . , Sub - 
aru/Hyper Suprime-Cam (HSC) (|Mivazaki et al.l I2OO6I '). 
SLSNe powered by the interaction can b e detected up to z ~ 
5 (e.g.. lTanaka et al.ll2012l : ICookell2008l) where as PISNe ca n 
be detected up to z ~ 2 (e.g.. iPan. Loeb. fc KaserJ 120121 '). 
Upcoming NIR surveys by, e.g., James Webb Space Tele- 
scope 

(JWSTfl 

Eucliqj, or Wide-field Imaging Surveyor for 
High-redshi ft (WISHJj, can approach to much higher red- 
shifts (e.g., iTa naka e" al.l |2012|: |Pan. K asen. fc Loeb| [20il : 
iHummel et al. 2011; S cannapieco et a l. 2005). Here, we 
show how the LC of SN 2006gy is observed if it appears 
at high redshifts. This can be used for the identification of 
candidate transients obtained in upcoming optical and NIR 
surveys. 

Figure [T3] is the i band LCs of the model D2 [w = 5) 
and Fl (w = 0) for several redshifts. As the model Fl has 
similar color to SN 2006gy at around the LC peak, the Fl 
model provides better estimate for the high-redshift SLSN 
LCs. Figure [16] is the LCs with the F277W filteiQ which 
is one of the wide filters planned for JWST. The central 
wevelength of the F277W filter is at 2.77 ^lm. JWST can 



^ http : //www . j wst . nasa . gov/ 1 
^ http://scl.esa.lnt/euclid 
^ http: //www. wlshmlssion. org/ en/index .html 

* ,http: //www. stscl . edu/ j wst /Instrument s/nlr cam/ Instrumentdesign/f liters/ 
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Figure 15. i band LCs of the models D2 {w = 5, left) and Fl (w = 0, right) at several redshifts. The horizontal line is the planed i 
band limiting magnitude of Subaru/HSC Ultra-Deep survey at one epoch (26.6 mag). 
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Figure 16. F277W band LCs of the models D2 {w = 5, left) and Fl {w = 0, right) at several redshifts. 



reach 30 mag and is abl e to detect SLSNe beyond z = 
10. Scc Tana ka et al.l (|2012l ) for the details of observational 
strategics and estimated number of detections of SLSNe with 
upcoming optical and NIR surveys. 



7 SUMMARY AND CONCLUSIONS 

We have shown that the interaction between SN ejecta and 
dense CSM is a viable mechanism to power SLSNe such 
as SN 2006gy. The interaction in the dense CSM accounts 
for the huge luminosity and the long duration of the SN 
2006gy LC. Shock breakout within the dense CSM is a key 
for the understanding of the interaction-powered SLSNe. 
Our canonical models have Mej — 20 Mq, Ecj — 10^^ erg, 
and McsM = 18 Mq {w = 5) or 15 Mq {w = 0) where 
the CSM is assumed to have a density profile of p oc r^™. 
The corresponding average mass-loss rate of the progenitor 
is about 0.4 Mq yr~^ if we assume that the dense CSM 
originates from a 100 km s~^ wind. Our steady mass-loss 
models (w = 2) fail to explain the SN 2006gy LC. No ^^Ni 
is required to explain the early LC of SN 2006gy. 

It is difficult to break the degeneracy among Mcj, Eaj, 



and McsM. One can obtain constraints on the progenitor 
of SN 2006gy based on the efficiency, as the conversion ef- 
ficiency of the SN kinetic energy to radiation becomes high 
when Moj is comparable to or less than A/qsm. The pro- 
genitor of SN 2006gy should be a very massive star because 
McsM ~ 18 Mq or 15 Mq As the conversion efficiency is 
~ 50 % at most and the radiation energy emitted by SN 
2006gy is more than 2 x 10^^ erg, Ecj should be larger than 
4 X 10^^ erg. 

We have also examined the effect of multidimensional 
instabilities in the dense cool shell on the model LCs. Such 
instabilities are expected to reduce the amount of kinetic 
energy converted to radiation. Our LC modeling is based on 
a one-dimensional radiation hydrodynamics code in which 
the multidimensional instabilities are implemented only in 
an approximate way. We have thus explore the effect qualita- 
tively. Further studies on the multidimensinal effect(s), per- 
haps using three-dimensional radiation hydrodynamics sim- 
ulations, are needed for better understanding of SNe pow- 
ered by the interaction. 

Finally, we have provided predictions from our model 
for high-redshift SLSNe. We have calculated the optical and 
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NIR LC evolution in the observer frame. The results can 
be used for the identification of SLSNe appeared at high 
redshifts in the future transient surveys. 

The existence of the very massive CSM close to the 
progenitor which is required to explain the LC of SN 2006gy 
challenges the current understanding of the stellar mass loss. 
The better understanding of SNe powered by the interaction 
will lead us to the better understanding of the mass loss 
mechanisms of massive stars. This is a critical key to reveal 
the evolution and the fates of massive stars which are a 
fundamental component in the Universe. 
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APPENDIX A: EFFECT OF THE SMEARING 
TERM ON LIGHT CURVES OF PULSATIONAL 
PAIR-INSTABILITY SUPERNOVAE 

We show the results of our additional investigations on 
the dependence of the smearing term Bq on the results 
of LC calculations. We show the results obtained from 
the pulsational-pair instabil i ty SN models presented by 
IWooslev. Blinnikov. fc Heeeil (2OO7I ) and show that the ef- 
fect of Bq on LCs is similar to our LC models. This means 
that the effect of Bq appeared in our models is not unique 
to our models. 

The models presented here are based on the same 
progenitor model (Mzams = HO Mq) shown in 
IWooslev. BUnnikov. fc Hegeij |20o3). The shell ejected 
by the first pulsation (24.5 Mq) is caught up by 
the materials ejected by th e next pulsation (5.1 Mq). 
IWooslev. Blinnikov. fc Hegeij l|2007l ) have artificially in- 
creased the kinetic energy of the materials ejected by the 
second pulsation to obtain the better fit to the LC of SN 
2006gy and their best LC model has four times as much ki- 
netic energy as that of the original model. The models shown 
here have more kinetic energy. The kinetic energy of the sec- 
ond pulse is increased nine times and the kinetic energy of 
the second pulse is 6.5 x 10^^ erg. The rest of the model is 
the sa me as those discussed in lWooslev. Blinnikov. fc HegeJ 
l|2007fl . 
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Figure Al. LCs of a pulsational-pair instability model with 
Bq = 1. The kinetic energy of the second pulse is increased nine 
times as large as the original model. Note that slightly higher 
absorption (A^ = 1.68 mag) is applied in the observational data 
points than the value adopted in the other part of this paper 
(Ajj = 1.3 mag). 



Figures IATIIA2I and IA3l show the results of the LC cal- 
culations with different Bq. Comparing the results to those 
shown in Figure lTT] we can see that the dependence of the re- 
sults on Bq is basically the same. This is because the source 
of the luminosity in both the models is kinetic energy of 
the material coming from inside and Bq directly affects the 
conversion efficiency of the kinetic energy to the radiation 
energy. For the accurate modeling of the LCs resulting from 
the conversion of the kinetic energy to radiation, we need 
more investigations on the multidimensional effects during 
the interaction. 
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